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1 Introduction 


Among the extensive literature on long-range dependence, relatively few publi¬ 
cations are devoted to cyclical long-memory processes or long-range dependent 
random fields. However, models with singularities at non-zero frequencies are 
of great importance in applications. For example, many time series show cycli¬ 
cal/seasonal evolutions. Singularities at non-zero frequencies produces peaks 
in the spectral density whose locations define periods of the cycles. A survey 
of some recent asymptotic results f or cyclical long-ran ge de pendent rando m 
processes and fields can be found in Ivanov et al. ( 2013 ) and lOlenkoi ( 20131) . 

In image analysis popular isotropic spatial processes with singularities of 
the spectral d ensity at non -zero frequencies are wave, J-Bessel, and Gegen¬ 
bauer models. lEsneio et al.l ( 2014 ) investigated probabilistic properties of spa¬ 
tial Gegenbauer models. A realization of the Gegenbauer random field on 
100 x 100 grid is shown in Figure [T| 
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Fig. 1: Simulated realization of the Gegenbauer random field. 


This article studies minimum contrast estimators (MCEs) of parameters of 
the Gegenbauer random fields. The MCE methodology has been widely applie d 
in different statistic al frameworks (see, for example. lAnh et, al.1 ( 20041 . [20071) : 
WeiLin et al.l ( 2012 )). One of the first works which used a MCE methodology 
for the para meter estim ation of spectra l den sities of stationary processes was 
the paper bv lTaniguclh (1987). Guvon ( 19951) introduced a class of MCEs for 
random fields. Anh et al.l (120041) derived consistency and asymptotic normality 
of a class of MCEs fo r stationary pro cesses within the class of fractional Riesz- 
Bessel motion (see iAnh et al. ( 1999 1). Resul ts bas ed on the second and third- 
order cumulant spectra were given bv I Anh et al. ( 200 7). They also provided 
asymptotic properties of second and third-order sample spectral functionals. 
These properties are of independent interest, since they can be applied to study 
the limiting properties of non para metri c esti mators of processes with short or 
long-range dependence. WeiLin et al. (1201 2) applied the minimum contrast 
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parameter estimation to approximate the drift parameter of the Ornstein- 
Uhlenbeck process, when the corresponding stochastic differential equation is 
driven by the fractional Brownian motion with a specific Hurst index. 

A burgeoning l iterature on s patio- t emporal estimat i on has emerged in re¬ 


cent decades (see IBeran et al.l (120091) . IChan and TsaH (l2012h 


( 2001|) . IGuo et al.1 ( 2009h . Li and McLeod ( 1986 1. Reisen et al. ( 2006h , among 


Giraitis et al 


others). One of the most popular estimation too ls applie d was the maximum 
likelihood estimation method (MLEl. lReisen et al.l (2006) addressed the prob¬ 
lem of parameter estimation of fractionally integrated processes with seasonal 
components. In order to estimate the fractional parameters, they propose 
several log-periodogram regression estimators with different bandwidths se¬ 
lected around a nd/or betwee n the seasonal frequencies. The same methodology 
was used bv lLi and McLeod ( 198fill for fractionally differenced autoregressive- 
moving average processes in the stationary time series context. Several con¬ 
tributions have a lso been made for MLE of long memory spatial processes 
(see, for example. Anh and Lunnevl (1995 1). For two-dimensional spatial data 
the paper by Basu and Reinsel ( 1993h introduced a spatial unilateral first- 


order autoregressive moving average (ARMA) model. To implement MLE 
they provided a proper treatme nt to bor der cell values with a substantial ef¬ 
fect in estimation of parameters. IBeran et al.l d2009lf addressed the problem of 
the least-squares estimation of autoregressive fractionally i ntegrate d moving- 
average (FARIMA) processes with long-memory. Cohen and Francos! ( 2002 1 
investigated asymptotic properties of least-squares estimators in regression 
models for two-dimensional random fields. Maximization of the Whittle like¬ 
lihood ha s been also cons i dered i n the recent li teratu r e on the MCE (see for 
examp le. Chan and Tsail (12012 1. Boissv et al. ( 2005 1. Leonenko and Sakhnol 


( 200611 1. Leonenko and Sakhno (200(1) gave a continuous version of the Whit 


tie contrast functional supplied with a specific weight function for the estima¬ 
tion of continuous-parameter stochastic p rocesses, de riving the consistency and 
asymptotic normality of such estimators. Guo et all ( 2009 1 demonstrated that 
the Whittle maximum likelihood estimator is consistent and asymptotically 
normal for stationary seasonal autoregressive fractionally integrated moving- 
average (SARFIMA) processes. 

Parameter estimation of stationary Gegenbaue r rando m pro cesses was con¬ 


sidered b y numerous authors, se e , for example, Cray et, al. ( 1989 1. Chund 


1 996ai lb|,! Woodwa rd ahl( 1998lf . lCollet and Fadilil ( 2006h , McElrov and Holanl 


( 201211 . Gray et al. (198911 used the generating function of the Gegenbauer 


polynomials to develop long memory Gegenbauer autoregressive moving-average 
(GARMA) models that generalize the FARIMA process. GARMA models were 
estimated by applying the MLE methodology. I Chun 2 (Il996alf also applied 
this methodology with slight modifications based on the conditional sum of 
squares method. IChuna (I1996bh extended these results to the tw o-parameter 
context within the GARMA process class. Woodward et al. (|l998h introduced 
a fc-factor extension of the GARMA model that allowed to associate the long- 
memory behavior with each one of the k Gegenbauer frequencies involved. 
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In this article we restrict our consideration to the estimation of long- 
range dependence parameters. It is motivated in part by cyclic processes, 
for which pole locations are known. Also in some applications the spectral 
density singularity location can be estimated in advance. Various methods, in¬ 
cluding semiparametric, wavelet, and pseudo-maximum likelihood techniques, 
of the estimation of a singu l arity location were d iscus sed by, for example. 


Giraitis et ah (2001), and [Ferrara and Guegan 


[Arteche and Robinson 

( 200lh . 

This paper introduces and studies the MCE of parameters of spatia l Gege n- 


bauer processes. Specifically, analogous of continuous-space results by Anh et al 


([20040 are formulated for random fields defined on integer grids. The consis¬ 
tency and asymptotic normality of the MCE are obtained using a spatial 
discrete version of the Ibragimov contrast function. The article develops a 
methodology to practically verify general theoretical assumptions for consis¬ 
tency and asymptotic normality of MCEs for specific models. The results pro¬ 
vide a rigorous platform to conduct model selection and statistical inference. 

The outline of the article is the following. In section 2, we start by intro¬ 
ducing the main notations of the paper. Some fundamental definitions and 
the main results of this article, Theorem 1 and 2, are given in §3. Section 4 
consists of the proofs of the main results. Section 5 presents simulation studies 
which support the theoretical findings. The Appendix provides auxiliary ma¬ 
terials that specify for our case the conditions that ensure the consistency and 
asymptotic n ormality of the M CE based on the Ibragimov contrast function 
formulated in Anh et al. ( 2004h . 

In what follows we use the symbol C to denote constants which are not 
important for our discussion. Moreover, the same symbol C may be used for 
different constants appearing in the same proof. 

All calculations in the article were performed using the software R ver¬ 
sion 3.0.2 and Maple 16, Maplesoft. 


2 Gegenbauer random fields 


This section intro duces some of the m ain defini tions of the Ge ge nbauer ran¬ 
dom f ields g iven in Espeio et al] (12014 1 (see also, ChuneJ ( 1996al lbh Gray et al. 
rilflSflh . and lWoodward et al . dl998h . :’or the temporal case). 

Let Y tl: t 2 , £ Z 2 , be a random field defined on the grid lattice Z 2 . 

Consider the fractional difference operator defined by 


V“ = (/ - 2 uB + B 2 ) d = (1 - 2 cos i *B + B 2 ) d = [(1 - e lv B){ 1 - e~ w B)\ 


(1) 


where B is the backward-shift operator, u = cos^, i.e. v = arccos(zt), |zt| < 1, 
and <1 £ (- 5 , 5 ) ■ Assume that Y satisfies the following state equation 

° = (/ - 2U1B, + Bl) dl o (/ - 2 u 2 B 2 + B 2 2 ) d2 Y tut2 = e tlM , 

( 2 ) 
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where V*, i = 1,2, is given by equation ©, with Bi , i = 1,2, denoting the 
backward-shift operator for each spatial coordinate, i.e. BiY tl ^ 2 = Y tl - i^ 2 , 
and B 2 Y tlt t 2 = Here, £ti,t 2 , (ti,t 2 ) G Z 2 , is a zero-mean white noise 

field with the common variance Me 2 + 1 = q~ 2 . The ra ndom field Y is called a 
spatial Gegenbauer white noise in lEspeio et al.l ((2014). 

By equation (© the Gegenbauer random field Y can be defined in terms of 
the inverse of the operator oVj 2 expanded in a Gegenbauer polynomial 
series as follows 


oo oo 






m=o U 2—0 


oo oo 


= E E 2 ) etl - niit2 _ na> 

ni— 0 ri2=0 

-y(d)/ 


(3) 


where di ^ 0, i = 1, 2, and Cn^ (it) is the Gegenbauer polynomial given by 

” ( } } fc!(n — 2 k)\r(d) ‘ 

The generating function for the Gegenbauer polynomials is given by 

OO 

Y J C^\u)b n = {l-2ub + b 2 y\ | 6 | < 1 , 

n— 0 

which explains the expansion of the inverse operator in © . 

In general, a random field Y is called invertible if the white noise St 1 ,t 2 , 
(ti,t 2 ) 6 Z 2 , can be expressed as the convergent sum 

OO OO 

£ti,t2 = ^ ^ ^ ^ ^ 1 , 712 ^* 1 — ni,t2— ri2 1 

ni=0 n 2 —0 

where E^oE^o l 6 n lf n a | < oo. 


^—* / i. i—u ^u 1 i 

_The Gegenbauer random field has the following property, see lEsneio et al 
(l2014h . 

Proposition 1 If 0 < di < \ and |it,| < 1, i = 1,2, then Y is a stationary 
invertible long range dependent random field. 

The spectr al densi ty of a stationary Ge genbauer random field is given by 


Ghune ( ~1996al b), and Hsu and Tsai ( 2009h . for the one-parameter case, and 


Esoeio et alJolldh . for the two-parameter case: 


/(A, 9) = Eh 1 1 _ 2 Ul e~ iXl + e~ 2iXl I 2dl 11 - 2u 2 e" iA2 + e“ 2iA2 1 


(27r) 


{|2cosAi - 2zti|} 2dl {|2cosA 2 - 2u 2 |} 

{2ny 


1 — 2^2 


(4) 
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where 9 = ( u,d) = (ui, u 2 , di, 0 ^ 2 ) £ 9 = (—1, l) 2 x (0,l/2) 2 , m = cos 
and — 7 r < A; < 7r, i = 1,2. Using the spectral density function © , one can 
compute the auto-covariance function of Y as follows: 

lUi,h,9) = ^pr(l-2dj)[2sin(^)]5” 2dl 

i— 1 


■/* 2 •/* 2 


w here P„{z) is the associated Leg endre function of the first kind, consult §8 


Abramowitz and Stegunl (|l£22 


From 


Chund 


Gray et al. (1 980h and Gradshtevn and R.vzhik 


( 1980ll . the following asymptotic approximation of the autocovariance func¬ 
tion can be obtained 


7(ii,i2,0) = n 




i= 1 


7 r sm 


2 di 


00 


sin(d i 7r)r(l - 2 di) c°s{j i i' i )^p-_ t 1 + X )]- 


The random field Y is long range dependent as its auto-covariance func¬ 
tion satisfies the condition j 2 )ez 2 l 7 (ii>j 2 >®)| = + 00 . A detailed discus¬ 
sion on relations between local specifications of spectral functions and the tail 
behaviour of au to-covariance functions of lon g range dependent random fields 
can be found in Leonenko and Olenkol ( 2013 1. 

Figure [2] gives an example of the spectral density and the auto-covariance 
function of the Gegenbauer random field for the values of the parameters 
u\ = 0.4, U 2 = 0.3, d\ = 0.2, and (I 2 = 0.3. 



Fig. 2: Spectral density and auto-covariance function for u = (0.4,0.3) and 
d= (0.2, 0.3) 


3 Asymptotic properties of MCEs 


Very detailed discussions about esti mation of parameters of seasonal/cyclical 


long memory time series were given in Arteche and Robinsonl (2000). Gir aitis et al 


(20011. iFerrara and Gueeranl (1200111 . and references therein. It was shown that 
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the case of spectral singularities outside the origin is much more difficult com¬ 
paring to the situation of a spectral pole at 0. In particular, estimators of 
parameters u,; and di have different rates of convergence and their joint distri¬ 
butions are still unknown under general conditions. 

Ferrara and Gueganl (2001!) pointed out that, in practice, parameter esti¬ 


mation of seasonal/cyclical long memory data is done in two steps. The first 
step consists in estimation of singularity locations, which does not represent 
a difficult issue. Then, the obtained values of location parameters are used in 
estimators of the long memory parameters. It has been shown that this 2-steps 
method provides quasi-similar results to simultaneous procedures. Similarly, 
the results of this paper may also be used in a hierarchical modeling frame¬ 
work. Namely, locations of the singularities can be included in the first level of 
the hierarchical procedure. Then, on the second step of the hierarchical pro¬ 
cedure, the long-range dependence parameters can be estimated by the MCE 
method, conditioning on the locations obtained on the first step. It also allows 
to construct a suitable weight function before applying the MCE methodol¬ 
ogy. Therefore, in this paper we only concentrate our attention on long-range 
dependence parameters. 

Suppose that the conditions imposed in Proposition [T| to ensure station¬ 
ary, invertibility and long-range dependence hold. Assume the value of the 
parameter u = (111,112) is known a priori or was estimated on step 1. Then 
9 = ( 01 , 62 ) = (di,d 2 ) G 0 = (0,1/2) 2 is the vector of parameters to estimate 
of the Gegenbauer random field defined by equation © (see also equation © 
for the corresponding spectral density). 

Let T tljt2 , ti,t 2 = 0be a part of a realization of the Gegenbauer 
random field. 

Let w( A), A e [—7r, 7r] 2 , be a nonnegative function. Suppose the condition 
A3 in the Appendix holds. We define 


a 2 ( 6 )= /(A, 9)w(\) dX 

J [—7T,7r] 2 


( 5 ) 


and consider the factorization 

f(\,9)=a 2 (9)V(\9). (6) 

For all 6 £ 0 the function >F(A, 6 ) has the following property 


[ <F(A, 9)w(X) d\ = 1. 

J [ —7T,7r] 2 


( 7 ) 


Let K( 6 o, 6 ) be a non-random real-valued function, usually referred as the 
contrast function, given by 

K(9 0 ,9) := / /(A, 9 0 )w( A) log d\, 

J [—7r,7r] 2 “'(A 
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and let the contrast field be 



/(A,floMA)log!?(A,0) dX, 


where 9q = (dio,d 2 o) is the true parameter value. In what follows, Pq denotes 
the probability distribution with the density function /(A, 9q). 

The empirical version Ut(9), T G Z, 9 € 0, is defined by 


U T (0) := — f I T {\)w(\)log&(\,6) dX, ( 8 ) 

J [ —7T,7r] 2 

where It(A) is the periodogram of the observations Yt ljt2 , ti, <2 = 0, of 
the Gegenbauer random field, that is, 


MA) := 


1 

(2ttT) 2 


^2 e ~ i(tiXi+t2X2)Y tut 2 
t± =0 £ 2=0 


The MCE is defined by the empirical contrast field Ut(9) and the contrast 
function K(9q,9) being K(9q,9) > 0, and having a unique minimum at 6 = 

Go- 

In particular, we choose 

w{X) = |A? - i/?| 01 |A| - v%\ a2 w 0 (X), X = (A,, A 2 ) G [— 7 r, 7 r] 2 , (9) 

where aj > 1, i = 1,2, u>o(A) is a positive function with continuous second 
order derivatives on [—tt,tt] 2 . 

Notice that by © the weight function w( A) is nonnegative and symmetric 
about (0,0). Due to the boundedness of w( A) the product w(X)f(X,9) is 
integrable. Thus, the choice of w( A) fulfils condition A3 in the Appendix. 
The developed methodology is readily adjustable to other classes of weight 
functions. 

Theorems |T] and [ 2 ] give consistency and asymptotic normality results for 
the MCE. 


Theorem 1 LetY tl j 2 , (£ 1 ,( 2 ) G Z 2 , be a stationary Gegenbauer random, field 
which spectral density satisfies equation (U). //C/t(0) is the empirical contrast 
field defined by equation m, then 

— y tl , t 2 satisfies the conditions A1-A6 in the Appendix; 

— the minimum contrast estimator 9t = (di,^) = arg mingle Ut{Q) € 0 
is a consistent estimator of the parameter vector 9. That is, there is a 
convergence in Pq probability: 

9 t A 6> 0 , T —*• 00; 

— dj, <t 2 (0o), T — y 00, where the variance estimator is given by 

d 2 T = f 7 r (A)w(A) dX. 

J [— 7T,7r] 2 






Rosa M. Espejo et al. 


Notice that Assumption A4 is not required to prove the last two statements 
in Theorem [I] but it will be used in the proof of Theorem [2 

To formulate Theorem [2] we introduce the following notations. The unbi¬ 
ased estimator of the correlation function 7 (^ 1 , t 2 ,0), £ = (£ 1 ,^ 2 ) £ Z 2 , of the 
Gegenbauer random field Y tlt2 is 


1 T—\t\\ T—\t2\ 

lT(t) — _ |^ 2 |) + 

Note that all indices of the random field in the sum above are within the set 
{(ii,i 2 ) : ti,t 2 = 0 , ...,T}, where the observations are available. 

The unbiased periodogram is given by 


T-l 


T-l 


JRAr.Aa) — 




E E « 

t 1 =l-Tt 2 = l-T 




} 7 r{t), 


and the corresponding empirical contrast field is 


U*(0) = - ( /*(A)u>(A)log!P(A,0) dX. 

J [ —7T,7r] 2 

We also define crfi = I^(X)w(X) dX and the associated adjusted MCE 


6 t = (^ 1 ,^ 2 ) = argrnin tj^{9). 


( 10 ) 


Theorem 2 If Y tl ;t2 , (ti, t 2 ) £ Z 2 , is a stationary Gegenbauer random field 
which spectral density satisfies equation with (di,d 2 ) £ (0,l/4) 2 , then 

— Y tllt , satisfies the conditions A1-A9 in the Appendix; 

— the adjusted MCE defined by (Hhd is asymptotically normal. That is, 


T(0 t - do) -A W 2 (0, S^(9o)A(9 q )S~\ 9 q )), T — > 00 , 


where the entries of the matrices S(0) = ( Sij{9 )) and A.{9) = ( aij(9 )) are 

J(0) = [ /(A, 9)w( A) -f— log l?(A, 0) dX = a\9) ( w(X) 

J\- TT.TTl 2 OVidtij J [ — . 7f 1 2 


a 2 1 

*(A,0)- 


d(9, : <9(9, 


<P(A,0) 5(9. 


9 <T(X,9)^-<E(X,9) 


5(9, 


dA, 


( 11 ) 


a y (5) =8 tt 2 [ / 2 (A,0) W 2 (A)Eiog(^(A,0))Eiog(^(A,0)) 

0 [ —7T,7rl 2 dVi 


dX 


= 8n 2 a 4 (0)J 


w 2 (X)^-T(X,0)^-T(X,0) dX. 

7r,7r] 2 OVi OUj 


( 12 ) 
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To avoid the edge effect I An 


1 et al 


odogram approach suggested by 


Guvod (198: 





employed the modified peri- 
We use their assumptions in 


Theorem [21 Note that some authors pointe d few problems in using , see 


IVidal-Sana l 20091) . lYao and BrockweU ( 2006h . and references therein. Various 
other modifications, for example, typed variograms, smoothed variograms, ker¬ 
nel estimators, to reduce the edge effect hav e been proposed. I t would be inter¬ 
esting to prove analogous of the results bv lAnh et al.1 (20041) and Theorem [2] 
for these modifications, too. However, it is beyond the scope of this paper. 
Moreover, it remains as an open problem whether the edge-effect modi f icati on 
is essential for the asymptotic normality or not, see lYao and Brockwell ( 2006 ). 


4 Proofs 


To prove the theorems we will use the following differentiability lemma. 


Lemma 1 Schilling (12005 . Theorem 11.5)] Let (X, F, p.) be a measurable 
space, A £ T, and 0 C R be an open set. Suppose the function F :Xx6>->-1 
satisfies the following conditions: 


1. For all 6 £ <9 : F(-,0) £ Ti(A); 

2. For almost all x £ A the derivative dF gg’^ exists for all 9 £ 0; 


3. There is an integrable function g : X —> R such that 
almost all x £ A. 

Then there exists J A F(x,9)dp(x) = f A 9F gg’ 0 ' 1 dp{x). 


dF(x,8) 

de 


< g{x) for 


Lemma 2 The function cr 2 (0) is bounded and separated from zero on 0. 
Moreover, its fii'st and second order derivatives are bounded on 0 and can 
be computed by 


d_ 

dQi 


w. 


[ — 7T,7r] 2 

=-2 r 


ixAfiXUulX 


' [ —7T,7r] 2 


log |2 cos A i — 2m\ w(X)f(X , 9) dX, 


(13) 


/* ^2 

^-^V(0) = w(X)— e 


-f(X,9)dX = A / log|2cosAi — 2m 

i J[—K, 7T] 2 


x log |2cosAj — 2uj\ w(X)f(X,0)dX , 


(14) 


where i,j = 1,2. 

Proof By the choice © of the weight function we obtain 

sup f(X,9)w(X) < Too 

[— 7 T, 7 r ] 2 X & 


( 15 ) 












































10 


Rosa M. Espejo et al. 


and 


supcr 2 (0) = / sup/(A, 9)w(\)dX < 4 tt 2 sup /(A, O)w(X) < + 00 . 

0(zO J[—-K,tt] 2 0(z© [—7r,7r] 2 x 0 


,2 d 2 


Hence, cr 2 (0) is bounded. 

Note that sup[_ 7r 7r ] 2 x6 i {|2cosAi — 2ui|} 2dl {|2cosA2 — 2u 2 |}^ a2 < +°o- 
Also, by the choice of the weight function, there exists S > 0 and a set A 0 C 
[—7T, 7r] 2 of non-zero Lebesgue measure such that w( A) > S for all A G A 0 . 
Therefore, 


inf cr 2 (0) = inf / /(A, 9)w(\)d\ 

^ S\(A 0 ) 


sup[_ W 7r] 2 xe {|2cosAi - 2ui|} 2dl {|2cosA 2 - 2u 2 |} 

which means that er 2 (0) is separated from zero on 0. 

Now, to study -A-a 2 (9) we compute ^-/(A, 6), i = 1,2: 


1 2(^2 


> 0 , 


8 


ae. 

t 2 


d 


^-/(A,0) = 7^2 {|2cosAj -2uj\} 2dj -^-{|2 cos A*-2^1} 2di 


2 

= log ^(2 cos Ai - 2ui) -2 ^ [2(cos A,; - )] 


1 - 2 ( 2 * 


x [2(cos Aj — Uj)] j = — 21og |2 cos A* — 2m\ /(A, 9), (16) 


where j ^ i and j = 1,2. 

Using © and m we conclude that 


sup 

- 7 T, 7 r ] 2 X © 


u>(A)^-/(A,0) 


< +oo. 


(17) 


Thus, by © and Lemma |T] there exists 

8 


88 


■a 2 (9) = 


f w(X)^f(X,9) dX 

J f—7r,7rl 2 at, i 


and 


sup 

e 




< 47r 2 sup 


»(A)A/(A. 9 ) 


< +oo. 


[— 7 T, 7 r ] 2 X0 

It is not difficult to find ae d g g . /(A, 0). By (flBl) . for i, j = 1,2, i ^ j, the 
second derivatives of / are given by 

o2 o 

^a/( A,0) = —2 log 12 cos Ai -2uj| —/(A.0) = 4 (log 12 cos A^ - 2ui|) 2 /(A, 0), 


a 2 


<96» ? <90 ?: 


(18) 

/(A, 0) = 4 log |2 cos Ai — 2rti| • log |2 cos Aj — 2uj| /(A, 0). (19) 
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It follows front (1181) . (1191) . and © that 

d 2 


sup 

[ — 7 T, 7 r ] 2 X 0 




< +oo, i,j = 1,2. 


( 20 ) 


Finally, by m and Lemma |T| there exists 

B 2 I" a 2 

urA<»=\ 

and 


!W = l,.r“’ (A) 4k /(A ' 9) ‘' A 


sup 

e 


d 2 


dOjddi 


a 2 {9) 


< 47r 2 sup 

[ — 7 T, 7 r ] 2 xO 


”< A) 4k /(A -' ) 


< +oo. 


( 21 ) 


Proof of Theorem[ T] We will prove that the conditions A1-A6 in the Appendix 


are sa tisfied. Therefore, we will be able to apply Theorem 3 by Anh et al 
(2004) and obtain the statement of Theorem [TJ 


The condition Al holds, since 6 0 belongs to the parameter space 0 = 
(0, |) which is an interior of the compact set [0, . 

It follows from representation (jTj of the spectral density that /(A,0i) ^ 
/(A,0 2 ), for 0 1 ^02- Thus, the condition A2 is satisfied. 

The class of non-negative weight functions w(X) defined by © consists of 
symmetric functions. Note that 


| cos(Ai) — cos(i'j)! = 2 


Sill 


A i + Vi 


sm 


A i 


C |A?-!/?|, 


when Xi —>• ±Vi. Thus, by © and representation © of the spectral density 
we get w(\)f(\,0) G Li([— 7T, 7 r] 2 ) for all 0. 

To verify A4, that is, to prove 


V e 


[ &(\,0)w(\) d\= [ Ve<F(A, 9)w(X) dX = 0, 

J f— 7 T, 7 rl 2 J\—n,n ] 2 


we find 


»(A)Aa>(A, S ) = ^[A/(A. S ) 


w( A) 


a 4 (9) [ d0, 


d 




/( X,9) ( 22 ) 


Si (A,61) 


S 2 (A,6I) 


and apply Lemma 1. 

The same symbol C is used for different nonessential constants appearing 
in the calculations below. 

By © and the choice of the weight function we obtain 

sup I log |2 cos Xi — 2m\ | /(A, 9)w(X) < +oo. (23) 

[ — 7 T, 7 r ] 2 X O 
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Therefore, by equation (IIGil and Lemma [21 


|Si(A,0)| < 


C 


< 


c 


e Ll ([-TT,Tr] 2 ). 


a 2 (0) min© cr 2 (0) 

Now, by m, m and Lemma [2] we can estimate S 2 (A, 0) as 


\s 2 (\,e)\ = 


w{ A) 


^a 2 (0) /( X,0) 


08 


a\8) 


— min<9<7 4 (0) — 


(24) 


(25) 


Finally, A4 follows from (IM1) . (Ell , and Lemma Q] with g{x ) = C. 

Note that Li([—7r,7r] 2 ) D L 2 ([—7r,7r] 2 ) = L 2 ([—7r,7r] 2 ). To verify the con¬ 
dition A5 for the weight function w(-) we have to show that /(A, do M A)' 
log*F(A, 0) £ L 2 ([—tt, 7t] 2 ), for all 0 £ 0. By ([3]) and © the product /(A, 0 O )- 
w(A) log<F(A, 0) is bounded for all A except {A : A* = ±z/j, i = 1,2}. Let 
di = max(d,, d,o). Then, for \ —> ±z/j, i = 1,2: 


2 

[/(A, Oo)w(X) log<F(A, 0)] 2 < C JJ |A» ± z^ 2 ®*- 4 ^ log |A* ± £ Li([-tt, tt] 2 ). 

i—l 


Therefore, combining the above results we conclude that A5 holds. 

To verify the condition A6 we use the following function v(X) = |A 2 — 
i/ 2 |^|A 2 — /? £ (0,1/2). Note that 


1^2 - ^ 1 “^ 


log/(A, 0) ~ C 


1^2 - ^ 1 “^ 


(di log | A? - 1/?1 + d 2 log | A? -i/i|) -t 0, 


when \ Thus, by the choice of v(-) and properties of <r(0) the function 

h( A,0) = f (A) logiF(A, 0) = | Af - || A| - (log/(A,0) - 21og<r(0)) 


is uniformly continuous on [—7r, 7r] 2 X 0. 
Also, it holds 


/(A,0o) 


w{ A) 


«(A) 


< Cv 1 (A) £ L 2 ([-7r,7r] 2 ). 


Si nce the conditio ns A1-A6 are satisfied Theorem 1 follows from Theorem 
3 in Anh et al. ( 2004 1. ■ 


Proof of Theorem 2. To prove the asymptotic normality of the MCE in The¬ 
orem [2] we will show that the conditions A7-A9 of the Appendix hold. 

We begin by proving the condition A7. First, to verify the twice differen¬ 
tiability of the function <F(A, 0) on 0 we formally compute the second-order 
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derivatives of : 


d 2 


dOjdQi 


d 

&fW) 

a 2 (0)- 


/(A,0) 

00, 

a 4 (0) 




a 2 


a\0) 


ddjddi 

0 2 

00 ^ 00 ^ 


/( \ 0 ) 


\0) 


\0) 


/(A,0)- 




55/( a ’ 9) 


1 


a8(e) 


8_ 

M 


/(A,0) 


^( 0 ) - 


_d_ 

M 


ct 2 ( 0 ) 


/(A,e) 


Note that in LemmaHwe proved that the derivatives wk° 2 ^ 

6), and gg d gg . /(A. 0) exist. Hence, by the above computations and 
Lemma [5] the function & is twice differentiable on 0. 

In addition, by estimates (fTTll , (EUl) , (12T1) , Lemma© and the above represen¬ 
tation for gg . gg . A, 9) the product w(X)f(X , 9 0 ) gg d gg log?F(A, 0) is bounded 
on [— 7r,7r] 2 x 6>. Hence, for i,j = 1,2, 0 G <9 : 

w (X)f(X,9o) de d de l°gff(A,0) G Li([-7r,7r] 2 ) nL 2 ([-7r,7r] 2 ). 

To prove part 2 of the condition A7, we first note that by © and (USD: 


■•**<*' fl > 


—2 log 12 cos A i — 2m | 


a 4 (0) 


By Lemma©the second term is bounded. Hence, it follows from m and (P251) 
that the product w(X)f(X, 9 o)^- log<F(A, 9) is bounded on [— 7 r, 7 r ] 2 x 0, that 
implies 


w(A)/(A,0 o )^-log!F(A,0) G i fc ([-7r,7r] 2 ), 


k> 1, i = 1,2, 0G0. 


To verify the condition A 8 we first check the positive definiteness of the 
matrices 5(0) and A(0). 

The entries of S(0) can be rewritten as 




( 0 ) 



& W (X ,9) 


d 2 

06 , 80 , 


\og'F w (X,9)dX, 


where ^ W (X,9) = f(X,9)w(X)/cr 2 (9). 

By 0), ©, and Lemma [2] the function 'F,,, (A, 9) is integrable, i.e. there is 
a constant C such that \P W (X,9)/C is a density. Hence, 5(0) = Ca 2 (9)l(0 ), 
where I is the Fisher information matrix of the random vector X with the den¬ 
sity W w { A, 9) = & W (X, 9)/ /, 2 'I J W {X. 9)dX. Therefore, 5(0) is non-negative 







































14 


Rosa M. Espejo et al. 


definite. Note that 1(9) = — (e , where Qi = ^-<P W (X,9). 

The random variables Q\ and Q 2 are not a.s. linearly related which implies 
positive definiteness of 5(0) . 

The entries of A(0) can be rewritten as 

/ o o 

w2 (\)—V(\, 9)—E( A, 9) dX = Ccr 4 (0) E (Q,Qj ), 

-tt.ttI 2 °9i OVj 


where Qi = -^-E(X, 9) and the random vector X has the density j - d\ ■ 

As (E (QiQj)) i ■_ 1 2 is a non-negative definite matrix, A(9) is non-negative 
definite too. Moreover, it is positive definite, because the random variables Q\ 
and Q 2 are not a.s. linearly related. 

Now we compute elements of the matrix S(9). By (fill) 


Sij(9) = a 2 {9) f w( A) 

J [—7T,7r] 2 

= o*(0) 


5 2 1 

*(A,0)- 


dOtdOj 


° *(\9)4rH\9) 


E(\ ,0) 09f v '' ! “'50, 


d\ 


' [—7T,7r] 2 

w( x ) 

f (A, 0)cr 6 (0) 
d 



’&/( A,0) 

a 2 (0)- 


1 

S' 

1 

b 


_ 0 _ 

dOi 


/(A,0) a\9)- 


dOi 


\9) 


/(A,0) 


09, 


f(X,9) a\9)- 


50, 


a 2 (0) /(A, 0)1 IdA 


L 


w(A) 

u 2 (0) 


W (A) 

a 2 (0) 

d 2 


d 2 


00 j 09i 


/(A, 0) 


a 2 (0) + 


50 


■u 2 (0) 


50 * 


/(A, 0) 


50,50,; 


a 2 (0) 


/(A, 0) — 


_5_ 

50" 


a 2 (0) 


_5 

50^ 


/(A. 0) 


w(X) ( 5 2 

(w/ m 

^(A) 

a 4 (0)/(A,0) 

5 


A/(A.S) » 2 (9) - 




50, 


■/(A,0) 


a 2 (0)- 


50, 


■°’ 2 (0) 


/(A, 0) 
/(A, 0) 


50 


'/(A, 0) 


*( 0 )- 


_5_ 

50,' 


! (0) 


/(A,0) )dX. 


By (|5|) , (|13|) 5 (|14p , and m we obtain 


AiW = 3 [ 

J [ —7T,7T 

w(A) 


^(A) 

u 2 (0) 


50 


■^(0) 


50 * 


/(A,0) 


dA 


' [—7T , 7r] 2 /(A, 0) 






dA = 
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-4 


-7T,7r] 2 


log |2 cos Xi — 2m\ log |2cos Xj — 2uj\ w(X)f(X,9)dX. 


By (O the elements of the matrix A(9) are 


i,A e ) = 87r 2 cr 4 (6») J 


= 8tt 2 


[ —7T,7r] 

d 


[— 7T,7r] 2 

2 (A) 


w 2 (X )^(X,9)^(X,0)dX 


dOj- 


w 


d_ 

dQi 


/(A,0) cr 2 (0) — 


<96» 


■a 2 ( 0 ) 


/(A,0) 


30, 


/(A,0) <j 2 {6) — 


80 


-a 2 (0) f(X,9) ) dX. 


Hence, by (flGl) we get a^j{9) = Si — S 2 (i,j) — S 2 (j,i) + S 3 , where 


Si = 32 tt 2 ct 4 (0) 


' [ — 7T,7r] 2 


log 12 cos A i — 2m \ log 12 cos A j — 2uj\ w 2 (X)f 2 (X , 0)dA, 


S 2 (i, j) = 16tt 2 ct 2 (0) 




' [ —7T,7r ] 2 


log|2cosAi — 2ui\ iu 2 (X)f 2 (X,0)dX , 


S 3 = 8 tt 2 




^ 2 <’> 


' [ —7T,7r] 2 


w 2 (X)f 2 (X,0)dX. 


The proof of the condition A9 is based on the approach in lBentkusI (Il972l l. 
Notice, that by © there exists a factorization w( A) = u>i(A) • wj 2 (A) of w( A) 
such that both /(A, 0 O ) = /(A, 0o)uh(A) and 1172 (A) are bounded functions 
of A. For example, one can select the function u>i(A) to be equal a product 
of |Ai — vi \ 2dl \\ 2 — V 2 \ 2d2 and a positive smooth function on [— n, 7r] 2 . Let us 
denote 102 (A, 0) = w 2 (X)-^- log*F(A, 0). Then, 

t[ (EI^X)-f(X,9 0 ))w(X)^-log*(X,9)dX = T f {EI* T ( A)uii(A) 

J[-ir,7r] 2 d[-7T,7i-] 2 

—/(A, 9q))w 2 {X, 9)dX = T [ (EI^{X) Wl {X) ~ EI^(X))w 2 (X,9) dX 


(El? (A) — /(A, 0 o ))iS 2 (A, 0) dA, 


(26) 


where /£(A) denotes the unbiased periodogram of the random field with the 
spectral density /(A, 0o)- 

Notice that f( X. 9u ) is bounded. Hence, by the first statement of Propo¬ 
sition 2 in IGuvoiJ dl982 ) the last integral in (BUI) is 0 (T~ 2 ) and the second 
term in ECU) vanishes when T —» 00 . Therefore, to prove the condition A9, it 
is enough to show that the first term in (EHl) vanishes too. 

Let y(t, 9q) denote the auto-covariance function of the random field with 
the spectr al de nsity f(X,9o). By multidimensional Parseval’s theorem, see 
Brychkov et al. ( 1992h . we get 


[ (EI*(X)wi(X)-EI^{X))w 2 (X,9)dX 

J [— 7T,7r] 2 
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1 

( 2 ^ 



E E 

ti — l—T t 2 —l—T 


T—l T-l 

- E E e- i(Xltl+X2t2) i(t,0 0 ))ti 2 (\,6)d\ 

ti=l — T t 2 =l—T 

(wi(A) ^2 e~ l(Xltl+X2t2) -/(t,9o)I[i-T,T-i^(t 1 ,t 2 ) 

( tl ,t 2)€% 2 



_ E e l{Xltl+X2t2) ^{t,e 0 )I [1 _ TtT _ 1]2 (t 1 ,t 2 2jw 2 (\,9)dX 



f(x, #o)^t-i(A 


x ) dx 


f ( x , 0q)^t-i (A — x) dx ) w 2 (A, 0) dA 


f(x + A, 0 o )^t-i(®)(wi(A) - wi(A + a:)) dx)w 2 (X , 0) dA, 

( 27 ) 

where = (2 E-i)F ( an( £/ 2 ) /2) ) ( EiESEO is the Fe j® r 

kernel, dn-T.T-iF^iifo) is the indicator function of the cube [1 — T, T — l] 2 . 
Let f(\,0 Q ) £ L 2 ([—7r,7r] 2 ), i.e. (di,d 2 ) £ (0,1/4) 2 . Then, one can re¬ 


fine the approach in Theorems 2.1 and 2.2 by Bentkus ( 1972h for the two- 
dimensional case with p = 2. Namely, let us split the inner integral in (l27l) 
into two parts: the first integral is over the region A a := {(x\,x 2 ) : \x\\ < 
T~ a , \x 2 \ < T~ a } and the second integral is over [7r, 7r] \ A a . 

For a £ (2 max(di, d 2 ), 1/2), by the choice of u>i(A) and 


f $ T -l(x)dx = ( -E- -f 

U ’ 1 7r(T- 1) J T -a 


sin((r- l)n/2) V 

sin(xi/2) y 


dan 


< 


1 


7T 2 (T-1) 2 



dan 


sin(a:i/2)) 


C 


y2(l-a) ’ 


T 


the both above integrals are bounded by Cet/T , where e-r —► 0 when T — > oo. 
It implies that first term in (ini) vanishes and completes the proof of the 
condition A9. ■ 
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Fig. 3: Boxplots of sampled values Fig. 4: Sample probabilities 
of Ot- Po(\0t ~ 0q\ < e)- 


5 SIMULATION STUDIES 

In this section we present some numerical results to confirm the theoretical 
findings. 

Figure [3] demonstrates a series of box plots to characterize the sample dis¬ 
tribution of MCEs of the parameters di, i = 1, 2, as a function of T. To com¬ 
pute it Monte Carlo simulations of the Gegenbauer field with 100 replications 
for each T = 10,20,30,40,50 were performed. For the parameters u\ = 0.4, 
U 2 = 0.3, d\ = 0 . 2 , g ?2 = 0.3, and er^ = 1 realizations of Y tl} t 2 were simulated 
using the truncated sum En 1= o En 2 =o i n ©• For example, a realization of 
the Gegenbauer random field on a 100 x 100 grid is shown in Figure [1] We 
set the parameter values of the weight function in © to a\ = ci 2 = 2 and 
wq(X) = 1. The periodogram Fr(A) was computed and the minimizing argu¬ 
ment 6 t of the functional Ut{0) was found numerically for each simulation. 
Figure [3] demonstrates that the sample distribution of Ot converges to 0 0 as 
T increases. The plot of the sample probabilities Pq(\0t — 0 0 \ < e) in Figure Q] 
also confirms convergence in probability of Ot to Oq. 

For each generated realization we also computed the value of a^ using 
It (A). Analogously to Figures [3l []] plots in Figures [5j [ 6 ] support convergence 
in probability of &t to g 2 (0q) when T increases. Notice that by © we get 
a 2 (0 0 ) ss 74.736 for the selected parameters. The larger values of £ in Figure[ 6 ] 
comparing to Figure [4] are due to the difference in the scales for the parameters 
(small values measured in decimals) and variances (large values measured in 
tens). 

A * 

To verify the result of Theorem [2] we used sample values 0 50 which min¬ 
imized the functional U^(0) for each simulation. To avoid possible negative 
values the modified periodogram I^( A) was truncated at zero by the R pro¬ 
gram. Bearing in mind the edge effect and modified periodogram’s correction, 
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10 20 30 40 50 

T 


Fig. 5: Boxplots of sampled values 
of (jfi. 



Fig. 6 : Sample probabilities 
Pod^T ~ cr 2 (^o)| < e). 



Normal Q-Q Plot of di 


Normal Q-Q Plot of 62 



Theoretical Quantiles 


Theoretical Quantiles 


Fig. 7: Boxplots of sampled values Fig. 8 : Normal Q-Q plots for each 
~ ~ * 
of 0 50 . component of 0 5O . 


Figures [7] [5] demonstrate that the results are close to the expected ones even 

for the relatively small T = 50. The normal Q-Q plot of each component of 0 5O 

in Figures [ 8 ] matches with the theoretical normal distribution. To test the bi- 

variate normality hypothesis about 0 5O we used the Shapiro-Wilk, energy, and 

kurtosis tests of multivariate normality from the R packages MVNORMTEST, 

energy, and ICS. In all the tests, p-values (0.9491, 0.4605, and 0.5314) con- 
- * 

firmed that 0 T asymptotically follows a bivariate normal distribution. Simu¬ 
lations for other values of the parameters were run, with similar results. 

Hence, we conclude that the MCEs are consistent estimators and the dis- 
tributions of 0 T converge to the bivariate normal law. Note that the simulation 
studies not only comply with the obtained results for (di, tfo) € (0,1/4) 2 , but 
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also indicate that the theoretical results may be extended to all possible values 
of (di,d 2 ) in (0,1/2) 2 . 

6 DIRECTIONS FOR FUTURE RESEARCH 


The esti mation methodology based on the unbiased periodog ram was intro¬ 
duced in 
paper by 


Guvonl (Il982 . 1995 ) see also Hevde and Gav ( 1998li . Recently, the 


Robinson and Sanzi (120061) and the references therein provided a de¬ 


tailed discussion on the topic. It studied mainly the difficulties arising in the 
application of the methodology in high dimensions. In particular, they investi¬ 
gated problems arising in relation to non-uniformly increasing domain asymp¬ 
totics associated with different expansion rates of the studied domain in each 
spatial direction. In this paper, we considered the case d = 2 and restricted out 
attention to the case of uniformly increasing domain asymptotics. The case of 
non-uniformly increasing domain asymptotics is left for future investigations. 

An extended version of the derived results can be obtained for more gen¬ 
eral formulations of the unbiased periodogram. In particular, different growing 
rates can be allowed for each spatial dimension in the definition of the sampling 
area. For example, one can consider the following generalize d versio n ofAhe 
two-d imensional unbiased periodogram (see, for example, iRobinson and Sanzi 
( 120001 ) ) 

32(T)-1 


Ig{^ 1, A 2 ) — 


1 


sicn-i 

(27t) 2 ' 

v ’ *1 = 1 —Si(T)t 2 = l— S2 (T) 


E 




7r(*), 


where the functions gi(T), i = 1,2, satisfy some suitable conditions (for exam¬ 
ple, gi(T) —>■ oo, T —> oo, and gi{T) < CT , C < 1, for i = 1, 2, and sufficiently 
large T). 

An important area for future explorations is to extend the results of lAnh et al.1 
( 20041) and simultaneously estimate the locations of singularities and long- 
range dependence parameters using the MCE methodology. A feasible way to 
approach this problem would be relaxing the L2-integrability assumptions in 
conditions A5-A8. 

It also would be interesting to extend the methodology bv iBentkus ( 1972 :) 
to prove the condition A9 for all (di,^) in (0,1/2) 2 . 

Note that our simulation results show that the proposed minimum con¬ 
trast estimation methodology works in the case of uniformly increasing domain 
asymptotics. 


Supplementary Materials The codes used for simulations in this paper are 
available from the site https: //googledrive . com/host/ 
OB7UxM8o_bnBxdG9zNU9MdHFOQUU/MCE7o20Gegenbauer7 0 20f ields/ 


Acknowledgements The authors are grateful for the referee’s careful reading of the paper 
and many detailed comments and suggestions, which helped to improve the paper. Also, this 







































20 


Rosa M. Espejo et al. 


work has been supported in part by projects MTM2012-32674 of the DGI, MEC, and P09- 
FQM-5052 of the Andalousian CICE, Spain. A.Olenko was partially supported by the 2013 
LTU research grant. 


Appendix 


The conditions for consistency and asymptotic normality of the MCE for 
parameter s of stationary fractional Riesz-Bessel type random fields given in 


Anh et al. (12004 1 are specified below for random fields on 1?. 


Al. Let Y tl j 2 , t = (£1 ,* 2 ) £ Z 2 , be a real-valued measurable stationary 
Gaussian random field with zero mean and a spectral density /(A, 0), where 
A = (Ai, A 2 ) £ [— 7 r, 7 r] 2 , 6 £ 0, and 0 is a compact set. Assume that 
60 £ int( 0 ), where Oq is the true value of the parameter vector 8 . 

A2. If 0i 7 ^ 02 then /(A, 0i) ^ /(A, 02) for almost all A £ [— 7 r, 7r] 2 with 
respect to the Lebesgue measure. 

A3. There exists a nonnegative function w( A), A £ [—7r, 7 t] 2 , such that 

1 . w( A) is symmetric about ( 0 , 0 ), i.e. w( A) = u>(— A); 

2. io(A)/(A,0) £ L 1 ([- 7 T, 7 r] 2 ) for all 0 £ 0. 

A4. The derivatives Vg !f (A, 0) exist and it is legitimate to differentiate under 
the integral sign in equation 0, i.e. 


if (A, 6)w( A) dX= f V e if (A, 0)w( A) dX = 0. 

J [ —7T,7r] 2 

A5. For all 0 £ 0 the function w(X), X £ [— 7 r, 7 r] 2 , satisfies 

/(A, 9 0 )w(X) log if (A, 0) £ Li ([- 7 r, 7 r] 2 ) C\L 2 ([— 7 T, tt] 2 ) . 

A 6 . There exists a function v(X), X £ [— 7 r, 7t] 2 , such that 

1 . the function h(X,Q) = v(X) log if (A, 0) is uniformly continuous on 
[—tt, tt] 2 x 0; 

2. /(A, 9q)w(X)/v(X) £ L \([ 7 r, 7 r] 2 ) n L 2 ([-tt, tt] 2 ). 

A7. The function if (A, 0) is twice differentiable on 0 and 

1. /(A, 6o)w(X)-gJ^g- log!f (A, 0) £ L 1 ([-n,-rr] 2 )C\L 2 {[- tt, tt] 2 ), for alH, j, 
and 0 £ 0 ; 

2. /(A, 9q)w(X)-^- log if (A, 0) £ Lfc([— 7 r, 7 r] 2 ), for all i, 0 £ 0, and /c > 1. 
A 8 . The matrices S(0) = (Sy(0)) and A(0) = ()) with the elements 

defined by (HD and Q are positive definite. 

A9. The spectral density /(A, 0), the weight function w{ A), and the function 
log !f (A, 0) are such that for all i and 0 £ 0 : 



' [ — 7T,7r] 2 


{EIj’{X)—f(X, 9q))w(X)-^- logf (A, 0) dX —> 0, as T —> 00 . 
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